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A third post-Newtonian (3 PN) equation of motion for an inspiralling binary consisting of two 
spherical compact stars with strong internal gravity is derived under harmonic coordinate condition 
using the strong field point particle limit. The equation of motion is complete in a sense that it is 
Lorentz invariant in the post-Newtonian perturbative sense, admits conserved energy of the orbital 
motion, and is unambiguous, that is, with no undetermined coefficient. In this paper, we show 
explicit expressions of the 3 PN equation of motion and an energy of the binary orbital motion in 
case of the circular orbit (neglecting the 2.5 PN radiation reaction effect) and in the center of the 
mass frame. It is argued that the 3 PN equation of motion we obtained is physically unambiguous. 
Full details will be reported elsewhere. 

PACS Number(s): 04.25.Nx,04.25.-g 

Renewed attention has been paid to a high order post-Newtonian equation of motion governing inspiralling compact 
binaries in the context of the efforts for direct detection of gravitational waves [1,2]. It is well-known that detectability 
of the gravitational waves emitted by the binaries and quality of measurements of astrophysical information (e.g. 
masses) depend on accuracy of theoretical knowledge of the waveforms [1], and hence partly of dynamics of the 
binaries. 

The 3 PN approximation has been a subject of much discussion because of its ambiguity reported originally in 
Jaranowski and Schafer [3]. In fact, the 3 PN ADM Hamiltonian in the ADM- type gauge obtained in [3] has two 
undetermined coefficients (akinetic and w s t a tic) and the 3 PN equation of motion in the harmonic gauge derived by 
Blanchet and Faye [4] has one coefficient A undetermined within their framework. Both groups have used Dirac 
delta distributions, which cause divergences in general relativity, to express the point particles and inevitably they 
have resorted to mathematical regularizations. Damour et al. [5] pointed out that the undetermined coefficients may 
arise due to unsatisfactory features of the regularizations they have used in [3,4]. Indeed, using the dimensional 
regularization, the work [5] have succeeded in determining both of the coefficients, namely, astatic = 0, which means 
A = —1987/3080 via relationship established in [6]. (^kinetic is related with Lorentz invariance and was fixed in [5,7]. 
Blanchet and Faye have developed a Lorentz invariant Hadamard Patie Finie regularization [8,9] and do not have any 
ambiguity other than A.) 

In gravitational wave data analysis, the reduction of predictability of the equation of motion due to the undetermined 
coefficient can become a problem. In fact, the 3.5 PN phase evolution equation and luminosity [10] unfortunately 
have four undetermined coefficients, one of which is A. 

Theoretically, a use of Dirac delta distributions and inevitable regularization should be verified in some manner. 
The perfect (physical) agreements among the results obtained by various authors with various methods [11-13] give 
a direct theoretical confirmation of the 2.5 PN result first derived by Damour and Deruelle [14]. It is important to 
achieve 3 PN iteration without introducing singular sources to derive unambiguous result and support the previous 3 
PN works which have used Dirac delta distributions. 

Based on our previous papers [12,15], we derive a 3 PN equation of motion for two spherical compact stars in 
harmonic gauge without introducing singular sources. Instead, we apply the strong field point particle limit [16] to 
deal with strong internal gravity of the stars. Our derivation is satisfactory in a sense that the equation admits 
conserved energy, is Lorentz invariant, and is unambiguous. In this paper, we shall show both of the 3 PN equation 
of motion and an associated 3 PN energy of the orbital motion in the center of mass frame and in the case of circular 
orbit. 

Below, we shall explain briefly our yet another derivation of a 3 PN equation of motion. Since this method is 
different from others, we mention some details specific to our method at the 3 PN order. After deriving an invariant 
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energy of the binary orbital motion, we shall compare it with that derived by Blanchet and Faye and fix the A 
parameter. Full explanation of our method including computational details will be reported in [17]. See also [12,15]. 

We write explicitly the post-Newtonian expansion parameter, e, which represents the smallness of the orbital 
velocities of the stars. The mass scales as e 2 from the post-Newtonian scaling. Then the strong field point particle 
limit [16] is achieved by setting the radius of the star to scale at the same rate as its mass, e 2 . The scalings of the 
mass and the radius enable us to incorporate in the post-Newtonian limit (e — > 0) a limit of a regular point particle 
whose internal gravity, ~ the mass over the radius, is strong irrelevantly to e. 

We derive an equation of motion via surface integrals of the gravitational energy momentum flux going through a 
sphere surrounding the star. For this method, we prepare two spheres Ba(t) = {x k \\x — za{t)\ < cRa} (A labels 
the two stars) on the r = constant surface, where r is the time coordinate in the near zone. Ba, called the body 
zone, is centered at the star A's representative point z a (t) and has a radius eR A where Ra is an arbitrary constant 
and smaller than the orbital separation but larger than the radii of the stars. We make the body zone radius shrink 
proportionally to e in the near zone coordinate (r, x 1 ) to ensure that the field on the body zone boundary due to the 
star is obtained by multipolc expansion when the e zero limit is taken. 

The Z-th multipole moments of the star A denoted by I a {t), including its mass, are defined as volume integrals 
over Ba of A TT = —g(T TT + t^) + X TTa/3 ' ,ap where g, T^ u and t^ L are the determinant of the metric g^, the matter 
stress energy tensor and the Landau-Lifshitz pseudo tensor. x^ ual3 ,afi arises since we use the wave operator of the 
flat spacetime instead of that of the curved spacetime when we solve the harmonically relaxed Einstein equation. 

I L A = e 2 ( d 3 a A A TT a L A , (1) 

Jb a (t) 

where we introduced the body zone coordinate a A = e~ 2 (x l — z a (t)) and multi-indices notation L = i\ ■ ■ ■ i\ (I > : 
integer). a A = a A ■ ■ ■ a A . The mass ttla = lfm^o P A with P T A = I A so defined would be the ADM mass if the 
companion star were absent and the body zone radius is taken to be infinite in the body zone coordinate. 

In deriving a 3 PN equation of motion, it is important in our formalism to notice that the body zone B A is a sphere 
in the frame where the star A orbits, but B A is not a sphere in the generalized Fermi frame [20] where the star A 
is at rest and the effect of the gravitational field due to the companion star is removed (as much as possible) except 
for, namely, the tidal effect. We define stars to be spherical in the generalized Fermi frame. Now, we define the 
"intrinsic" multipole moments I a (t) on f = constant surface in the generalized Fermi coordinate (f,x 4 ) as a volume 
integral over a sphere B a (t) centered at the same world event and with the same radius cRa as Ba(t). For example, 
the difference between the symmetric tracefree quadrupole moments is defined as e 8 5I A l:,> = e 8 I A l:)> — e 8 I A l3> = 
j' BA ^d?y A h TT y A l y^ > — J Ba ^ d 3 y A A rT y A l y A > where < .. > denotes symmetric tracefree operation on the indices 
between the brackets. For a "spherically symmetric" compact stars, the "intrinsic" tracefree quadrupole moment 
J^' 3> vanishes. Up to the 3 PN iteration of the gravitational field, SI A t:>> arise mainly due to the Lorentz contraction 
and can be evaluated as a surface integral over dB A . 

51< ij> = e- 6 \v k A v l A I dSk&vf&hr + 0(6 4 ) = -e 2 |m»i> + 0(e 4 ). (2) 

A JdB A 

In our formalism, it is possible to derive the 3 PN field for an isolated star (by taking limit where the mass of 
the companion star is zero). SI A %:i> is a necessary term to obtain the correct expression of the 3 PN field for such a 
system. Although other multipole moments defined over B A possibly hide purely monopole terms, only the quadrupole 
moment is found to be relevant up to the 3 PN order. Clean extraction of monopole terms from the multipole moments 
defined in our previous works is a problem at the 3 PN order specific to our formalism. Blanchet and Faye on the 
other hand elaborate their generalized Hadamard Partie Finie regularization [8] in a Lorentz invariant manner [9], 
and have properly taken into account of special relativistic kinematic effects including the Lorentz contraction. 

Now, let us briefly explain our derivation of the equation of motion. The local conservation law of the total energy 
gives an evolution equation for four momentum of the star and relationships among the multipole moments, namely, 
momentum- velocity relation. The last read as P\ — P A v A + Q l A + e 2 dD A /dr where P A = e 2 J Ba d 3 a A A Tt and P A are 
the three momentum and the energy of the star A. Q\ = e~ 4 § Ba dSk{A Tk — v A A TT )y A arises since the (pseudo-)stress 

energy momentum of the field extends outside of the star [12]. Q l A can be evaluated explicitly and do contribute the 3 
PN velocity momentum relation. We can define the representative point z a (t) of the star A by specifying the dipolc 
moment of the star, D A = I A (, e.g., z a (t) corresponding to a condition D l A — may be called the center of mass of 
the star A from an analogy of the Newtonian dynamics). The relationship between the energy P A and the mass m A 
can be obtained by integrating functionally the evolution equation of P A , which is expressed as surface integrals and 
can be evaluated explicitly up to the 3 PN order, in the form as P A — m A [l + 0(e 2 )]. 
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Combining the mass energy relation, the momentum velocity relation, and the evolution equation for the four 
momentum, we obtain the general form of the equation of motion [12]; 



« T JdB A JdB A 

+e- A v\ (I dS k A kT -v\l dS k A 7 

\JdB A JdB A 



dQ\ 2 d 2 D\ . nT .di 



_^_A +(mA _ PI )_^ (3 ) 



Note that A^ u = (—g)t1 u L + y^ VOL ^ \ a p on 8Ba, since 8Ba is well outside the star by construction of the body zone. 
The acceleration in the right hand side of Eq. (3) should be understood to be lower order acceleration than in the 
left hand side. The terms in the right hand side of Eq. (3) are completely expressed as surface integrals over the 
body zone boundary except for D\ to be specified. The surface integral approach enables us to derive an equation 
of motion irrelevant to the internal structure of the star (Effects of the star's internal structure on the orbital motion 
such as tidally induced multipole moments appear through the field and hence the integrand A M "). The scaling of the 
body zone radius cRa ensures that we have an equation of motion for compact stars. 

The field equation coupled to the matter equations mentioned above is the integrated relaxed Einstein equation 
under harmonic gauge (W = rf v — sj— gg^" where rf v = diag(— e 2 , 1, 1, 1). The harmonic condition is then h^ v \ v = 0) 

Jc(t,x*) \x-y\ 

We split the flat light cone C(r, x l ) into four parts; two body zones Ba, near zone outside the body zones N/B 
surrounding the binary and the far zone outside the near zone. For Ba and N/B contributions to the field W v {t, x l ), 
we expand the retarded field about the near zone time r. Then multipole expansion of the star is used to evaluate 
the two body zone contributions. The N/B contribution is basically evaluated with help of super-potentials (a super- 
potential here means a particular solution valid in N/B of a Poisson equation with a non-compact support source in 
N/B). Unfortunately, it was not possible to find all the necessary super-potentials explicitly at the 3 PN order. For 
integrands in Eq. (4) such that we could not find super-potentials in closed forms, after making the retarded expansion 
we leave the Poisson integrals unevaluated and substitute the (not-integrated) field into A^ v in Eq. (3). Then we 
perform the surface integrals in Eq. (3) (with respect to the spatial variable x % in Eq. (4)) first and next perform the 
remaining volume integrals (with respect to the spatial variable y l in Eq. (4)). In other words, we extract the parts of 
the field necessary to derive the equation of motion by interchanging the order of the integrations in Eq. (3) and Eq. 
(4) . As a check, we applied this method on the integrands for which the necessary super-potentials can be derived in 
closed forms, and found that both methods give the same result. Finally, we have dealt with the integration over the 
far zone using the DIRE method [13,18] and found that it does not contribute to the 3 PN equation of motion [13,19]. 

Using the method mentioned above, we obtain a 3 PN equation of motion for a two spherical compact stars binary. 
We here present the 3 PN relative acceleration in the case of the circular orbit and in the center of mass frame, which 
is an appropriate equation to inspiralling binaries. 

dV i 

— = -Clt n r\ 2 + e 5 2 . 5PN A\ (5) 

where V 1 — v\ — v\ is the relative velocity and 2.5PN^ is the relative acceleration at the 2.5 PN order (the radiation 
reaction term). The 3 PN orbital angular frequency Qi D is, (for comparison, we adopt similar notations as in [4]) 

l + e 2 7 (-3 + ^)+ 6 4 7 2 (e+^u + u 2 
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m 2 2 n = 7 

a , / „ f 2375 41tt 2 „„, /r 12 \] 19 , 



+ 0(6 7 ), (6) 



where m = mi + m,2, v = rnxm-i/rn 1 , 7 = m/r\2, and lni?o = (mi ln(ei?i) + m 2 ln(ei?2))/m. In Eq. (5) with (6), 
the representative point of the star A, z\(t), is defined by setting D\ = c 4 5 A q = — 86e 4 m^iva^/9, where Nd A ls the 
Newtonian acceleration. This choice makes the three momentum P l A parallel to the velocity v\. We note that there 
is no arbitrary parameter other than the body zone radii ei?A in the 3 PN relative acceleration. We here note that it 
is not allowed to fix the A parameter by comparing Eq. (6) with the corresponding result of Blanchct and Faye [4], 
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since the harmonic condition both groups have used does not fix the gauge completely [17] and the expression of the 
3 PN orbital angular frequency in terms of the coordinate distance m/r^, Eq. (6), is gauge dependent. (From the 
same reason, we can not fix A using Eq. (13) below.) 

We can remove away cRa dependence from the 3 PN relative acceleration, Eq. (5), physically by a suitable 
redefinition of the representative points of the stars. In fact, by setting 



D 



A. Now 



e 4 S 



AQ 



e4 y m Aiv<4 In 



ri2 



eR A J 



(7) 



we obtain the 3 PN relative acceleration free from any arbitrary parameter 



2.5PN-' 



(8) 



with m 2 n 2 = m 2 fl 2 n - 22e b 7 z/ln(ri2/-Ro)- This observation in fact is the case in general cases ( i.e., in general orbits 
not in the center of mass frame); The 3 PN equation of motion in general cases we have derived is physically free from 
any ambiguity. 

A reason why we are concerned with In cRa dependence is the following. Blanchet and Faye have introduced four 
arbitrary parameters in their regularization procedure, two of which appear in the regularization of the field having 
two singular points, and the others appear in the regularization of equations of motion for those two points. They 
showed that the two of those parameters can be gauged away, while the other two were consumed to make their 
equations of motion conservative (modulo the radiation reaction effect), and they found there remained one and only 
one parameter, A, although relationship between energy conservation and regularization parameters associated with 
point particle description is not clear. Our redefinition of the representative points (7) corresponds to their gauge 
transformation. Then, their observation makes us check if it is physically allowed to remove the In cRa dependence 
in our 3 PN equation of motion, since we introduced only two arbitrary parameters cRa and we have no freedom to 
make our equation motion conservative by adjusting these two parameters if we remove them away. Thus, we have 
two problems to be solved in our method; removal of In eRA and an energy conservation. For lack of space, here we 
show some facts which support naturality of Eq. (7). The energy conservation problem will be addressed in [17], 
there we shall show our equation of motion and an associated conserved energy of the binary orbital motion in general 
cases. 

Let us consider the harmonic condition. 



= /i r % = 4e 4 



= /i v 



E 

A=l,2 

E 

A=l,2 



1 dP\ r% 



t a dr 



dSi 



dB A 



y- 1 dP\ 

A =i,2 rA dT 



P T A v\ 



f 2dD\ 



(A Tl - v A K TT ) + 



E 

A=l,2 



JdB A \x-y\\ A I 



+ ■ 



(9) 



(10) 



where 



irrelevant terms. These equations are manifestation of the fact that the harmonic condition is 



consistent with the evolution equation of P A , the momentum- velocity relation, and the equation of motion (and 

• "). Thus, if logarithmic dependence of cRa arises from the 



relations among higher multipole moments, hidden in ' 
equation of motion (essentially the second term of Eq. 



(10)), P\ must have the same logarithmic dependence (times 



minus sign) to ensure harmonicity. This and the momentum velocity relation in turn mean P A , v A — dz A /dr or D A 
have corresponding logarithmic dependence. We found that P A have no logarithm up to the 3 PN order. Therefore 
z\ or D A should have logarithms. This is consistent with the fact that a choice of D A determines z A . z A depends on 
logarithms if the old choice is taken, while it does not if our new choice is taken. 

The second fact which supports our interpretation is as follows. We find that the near zone dipole moment D l N 
defined by a volume integral of A TT y t becomes 



e 2 D l 



N 



■ [ d 3 yA T V = E P aA + e 2 Y. D A + ^ [ d 3 yA^y\ 

Jn A=l,2 A=l,2 Jn / B 



(11) 



Then if we take the old choice of D A , the volume integral becomes 
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d 3 yA TT f = e 



N/B 



4 22 V 3 i i ( r 12 

m^Ara A In 1 



,4=1,2 



+ 



(12) 



where terms denoted by 



are independent of Notice that the near zone dipole moment can be freely 



determined, say, D l N — 0, since we can define the origin of near zone freely in general [21]. By taking temporal 
derivatives of D l N twice, we see that D^ Now gives the definition of z\{t) in terms of which the 3 PN equation of 
motion is independent of cRa- 

Finally, we show the 3 PN conserved energy (neglecting the 2.5 PN radiation reaction force) of the circular orbital 
motion in the center of mass frame. Using Eq. (5), we have 



^In (7) 



+ e V I —FT + + T7 V 



l + e 2 7 

235 
"64 



-3 + l)+<V 



27 
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64 



49 
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10141 123tt 2 



576 



64 



3 \R 



+ 0{e T ) 



(13) 



In terms of x = (mQi n ) 2 > 3 we obtain the 3 PN energy in an invariant form 



E ln (x) = - 



mvx 



12 



u\x + e 1 



27 19 



24 



V ) x 2 



675 f 34445 205tt 2 



64 



+ 



576 



96 



] 155 2 

> v v A - 

J 96 



35 
5184' 



i/ 3 x s 



+ 0(e 7 ). 



(14) 



Similarly, using Eq. (8), we have E(j) = £q n (7) + (ll/3)e 6 r?w 2 7 4 hi(ri2/-Ro)- Here we note that the difference between 
E(pj) and E\ n ("f) is merely due to the redefinition of the dipole moments (or cquivalently, a coordinate transformation 
under the harmonic coordinate condition). The invariant energy E(x) is the same as E\ n (x) but with x = (mi!i n ) 2 / 3 
replaced with x = (mfl) 2 ' 3 . This third fact that the energy has the same form for both definitions of the representative 
points of the stars when we write the energy in terms of the orbital angular frequency which is an observable supports 
that the apparent body zone radii dependence of the 3 PN relative acceleration has no physical effect on the orbital 
motion. 

We have thus derived a 3 PN equation of motion which takes account of strong internal gravity and avoids any 
ambiguity. Comparing our result, Eq. (14), with the corresponding result in [4], we determine the coefficient unde- 
termined in the Blanchet and Faye 3 PN equation of motion as A = —1987/3080. This value of A is consistent with 
the result of Damour, Jaranowski, and Schafer [5]. Thus, our result (indirectly) validates their use of the dimensional 
regularization in the ADM Hamiltonian approach in the ADMTT gauge. Finally, we note that Blanchet et al. [22] 
have recently obtained the same value of A, who computed a 3 PN equation of motion in the harmonic gauge using 
dimensional regularization. 
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